Mon. Not. R. Astron. Soc. 000, [THT2] (2012) Printed 12 March 2013 (MN WT^ style file v2.2) 



Two new ultra-cool benchmark systems from 
WISE+2MASS 



'cT)' J. I. Gomes, ^ D. J. Pinfield/ F. Marocco,^ A.C. Day- Jones, ^'^ B. Burningham,^ 
^ Z. H. Zhang, ^ H.R.A. Jones, ^ L. van Spaandonk,^ D. Weights^ 



0^ 

o 

6 



> 

00 

(N 
(N 

rn 
O 
cn 



X 
S3 



Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Hertfordshire, ALIO 9AB. 
^ Departamento de Astronomia, Universidad de Chile, Camino del Observatorio 1515, Santiago, Chile. 



ABSTRACT 

We have used the 2MASS all-sky survey and the WISE to look for ultracool dwarfs 
that are part of multiple systems containing main sequence stars. We cross-matched 
L dwarf candidates from the surveys with Hipparcos and Gliese stars, finding two new 
systems. The first system, G255-34AB is an L2 dwarf companion to a K8 star, at a 
distance of 36 pc. We estimate its bolometric luminosity as log L/Lq = —3.78 ± 0.045 
and Toff = 2080 ± 260 K. The second system, GJ499ABC, is a triple, with an L5 dwarf 
as a companion to a binary with an M4 and K5 star. These two new systems bring the 
number of L dwarf plus main sequence star multiple systems to twenty four, which we 
discuss. We consider the binary fraction for L dwarfs and main sequence stars, and 
further assess possible unresolved multiplicity within the full companion sample. This 
analysis shows that some of the L dwarfs in this sample might actually be unresolved 
binaries themselves, since their Mj appears to be brighter than the expected for their 
spectral types. 
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1 INTRODUCTION 

In the last years an increasing number of binary systems 
with ultracool components have been discovered and stud- 
ied. Ultracool dwarfs (UCDs) are defined as dwarf stars with 
spectral types of M7 or late r, and inc lude substellar brown 
dwarfs ()Baraffe et al. 1998i ; |Jones fc Steele 2001j : masses < 
0.072 M0) and some very low mass stars. UCDs populate 
the lower temperature rauRe, from effecti ve temperatures 
Teff ~ 2400 K to 500 K (|Kirkpatrick IbOOST ). 

Beyond the late M dwarfs, lies the L dwarf se- 
quence, which have Teg ranging from 2250 K to 1400 K 
l|Kirkpatrick et cil~lll999l ). While for earlier types there can 
be some hydrogen burning very low mass stars, most L 
dwarfs will be substellar-mass brown dwarfs. The sequence 
ex tends into the coole r T dwarf regime and to T^s of ~ 500 
K jLucas et al. II2OI0IV The Y dwa rf spectral type, first pro- 
posed bv lKirkpatrick et alTI (Il999l') has only recent ly seen its 
first confirmed objects with lCushing et al. I (|201ll ) discover- 
ing six Y dwarfs from the Wide-field Infrared Survey Ex- 
plorer (WISE). More recently lKirkpatrick et al. I l|2012l ) pre- 
sented seven more Y dwarfs from WISE, bring ing the current 
total of confirmed objects to 13. The WISE (|Wright et al. I 
I2OIO I) satellite has surveyed the entire sky in four infrared 
bands ( Wl, W2, W3 and Wi), at wavelengths of 3.4, 4.6, 12 
and 22 /im. Bands Wl and W2 give the best sensitivity to 



the coolest brown dwarfs due to the strong absorption cen- 
tred at 3.3 /im ( Wl) owing to methane for objects with T^s 
< 1500 K. At 4.6 nm {W2), however, there is no methane; 
as such the Wl — W2 colour is very sensitive to cool brown 
dwarfs. 

Several other large-scale optical and near-infrared 
(NIR) surveys have already proved to be effective tools at 
finding and studying large popu lations of UCDs: the Two- 
Micron All-Sky Survey (2MASS: ISkrutskie et al~ll2006l ) and 
the Deep Near Infrared Survey of the Southern Sky (DE- 
NIS; Epchtoin ct al. 19 9% in the NIR, ari d the Sloan Dig- 
ital Sky Survey (SDSS: lYork et al. I modi ) in the optical. 
As brown dwarfs cool and fade significantly over time this 
leads to a mass-age degeneracy, with lower mass younger 
brown dwarfs having similar Tcff and luminosity to higher 
mass, older ones. Measuring atmospheric properties (Toff, 
log g, [m/H]) of brown dwarfs in order to infer mass and 
age, is thus crucial. Using near-infrared and optical spectra 
and model fitting techniques, one can attempt to constrain 
physical properties for UCDs. However dust condensation, 
non-equilibrium chemistry and complex molecular opacities 
make it very challenging to accurately model such atmo- 
spheres, and the reliable fitting of properties with sp ectral 
models is not currently reliable l|Pinfield et al. |[2012l ). 

One way to overcome these problems is to identify 
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UCDs whose properties can be estimated in a relatively 
independent way. We refer to such objects as benchmark 
UCDs. They can be used as a testbed for prevailing theories 
and models. UCDs that are wide companions to main se- 
quence (MS) stars of the Galactic disc are fairly numerous, 
cover the full range of age and me tallicity, and are there- 
fore particularly useful benchmarks. Ivan Biesbroeckl (119611 ) 
was one of the first to search for these systems and the 
method used has bee n followed by many other successfull 
studies (e.g . Ixirkpat rick et al. I 2001: W ilson et al~l l200ll: 



Table 1 . Uncatalogued regions that have been removed from the 
initial sample 



Gizis et al. 



20011: iPinfield et al 



2001; W ilson et al. 
2006: B urningham et alTI 



20091: iFahertv et al I |2009| . I201Q : .Dav- Jones et a. I l201ll : 



Pinfield et al. II2012I ). 



Since the peak of the spectral energy distribution (SED) 
of L dwarfs occurs at NIR wavelengths, 2MASS and WISE 
are ideal to discover large numbers of these objects. In fact 
both surveys cover the whole sky and are well matched which 
makes them a powerful tool to detect L dwarfs. Despite their 
relatively large pixel sizes of 1 arcsec in 2MASS and 1.375 
arcsec in WISE, the 10 year baseline between these sur- 
veys also makes it possible to accurately determine proper 
motions. This is important as it allows confirmation of the 
binarity in these systems. 

We present here the measurement of proper motion for 
10 L dwarf candidates and the discovery of two new wide 
binary systems. Section 2 describes the technique used to 
identify these systems and how we assessed their common 
proper motion, describing our follow-up spectroscopic obser- 
vations whilst in section 3 we present the two new binary 
systems in more detail, addressing the issue of possible con- 
tamination. In section 4 we discuss results, estimate a binary 
fraction for these type of objects and investigate the possi- 
bility that some known UCDs with MS stars companions 
are actually unresolved objects. Finally our conclusions are 
presented in section 5. 



2 CANDIDATES SELECTION 
2.1 L dwarfs selection 

The UCD candidates were selected using NIR photome- 
try available from the 2MASS all-sky point catalogue. This 
database is ideal to search for L dwarfs as it covers the whole 
sky and the NIR JHKs filters sample the peak of the SED of 
such objects. Colour selection criteria were applied to obtain 
a first sample of L dwarf candidates. This selection required 
the following: 

0.5 ^ J-_ff 1.6 

1.1 ^ J - a: 2.8 

OA ^ H - K ^ 1.1 
J-H ^ 1.75 {H -K) + 0.37 
J -H ^ 1.65 {H - K) - 0.35 

We only consider candidates with a J band magni- 
tude ^ 16.0. These constraints were chosen based on the 
colours of known L dwarfs from Dwarf Archive^ To avoid 
source confusion and contamination from MS and giant stars 
some high density regions were excluded, such as the galac- 
tic plane (b ^ |25°|) and the Large and Small Magellanic 
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^ See |http: / /Dwarf Archives. org| 



Clouds. Two other reddened and overcrowded regions were 
also removed. Information about these regions can be seen in 
Table [1] The equatorial poles were also avoided {& |86°|) 
as 2MASS is not complete in these regions for optical cross- 
matching. The total area of the sky included in the first 
sample was 22, 178 deg^, or 53 per cent of the sky. 

We have kept objects with no optical counterpart within 
5 arcsec in the USNO-A2.0 catalogue or, if they do have 
counterparts, with R — K > 5.5. By doing this we ensure 
that all detections with optical counterparts are of spectral 
types later than M6 as these objects have extremely red 
colours in both the optical and IR domain due to their low 
Tcff. A signal-to-noise ratio (S/N) criteria was imposed to 
exclude low-quality photometric data (S/N > 5 for all the 
bands). Flags including ccflg=000, prox>6, gaLcontam=0 
and mpflg=0, were applied to make sure that the sources 
were unaffected by known artefacts, such as diffraction 
spikes from nearby bright stars, and to exclude extended 
sources and known minor planets. This first sample of pos- 
sible UCD candidates encompasses 28, 023 objects. 

We find that 269 of the 602 known L dwarfs were 
present in our initial list. The remaiiung ones did not make 
it through our initial cuts, either because they had J band 
magnitudes larger than 16, galactic latitudes outside our 
limit or they did not pass the quality fiags imposed. It is 
important to note that only 36 were excluded based on their 
colours, thus, the vast majority of L dwarfs have NIR colours 
inside the colour space we initially defined, thus validating 
our method. 



2.2 Selection of binaries 

2.2.1 Assessing photometry and separation 

In order to search for companions to our L dwarf can- 
didates, we selected a sample of M S stars from differ - 
ent catalog ues, namely the Hipparc os (|van Leeuwen IlioOTi ) 
and Gliese jGliese fc Jahreiss1ll99ll ) catalogues. We selected 
only Hipparcos stars within 50 pc as we only want to con- 
sider pairs up to this distance, whilst Gliese stars have a 
maximum distance of 25 pc (due to limitations of the cata- 
logue itself). The sample consists only of F, G, and K stars. 
Assuming the star's distance, we looked for L dwarf com- 
panions up to an on-sky separation corresponding to 20, 
000 AU, and angular separations of up to 10 arcmin, finding 
572 possible pairs. 

Next, we plot the UCD candidates on a colour- 
magnitude diagram (CMD), in order to assess if the their 
photometry was consistent with a real UCD at the com- 
panion's distance. Using the distance of the primary star. 
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Figure 1. Mj versus J — K colour-magnitude diagrams. The 
UCDs from the candidate pairs are shown as crosses. Inside the se- 
lection box, the filled triangles are the three known L dwarfs that 
turn up in our sample as possible binaries, the filled diamonds are 
UCDs candidates not confirmed and not in the dwarfarchive.org. 
The open triangles represent other UCDs with known parallax 
measurements from Dwarf Archives and are displayed for com- 
parison purposes only. 



we calculated Mj and with the J ~ K colour, obtained 
the plot shown in Fi g. [T] Following the method used 
by iPinfield et al. we selected pairs where the L 

dwarf candidate was within a specific region. This re- 
gion was defined using a sample of L dwarfs with known 
parallax, or with distances inferred from com panion stars 
l|Kirkpatrick et al. II2OOII: IWilson et al. Il2001^ and the Mj 
range from iKnapp et al.H ( 20041 ). We considered all of our 
pairs and using the CMD, identified 25 possible MS star- 
UCD binary systems. 



2.2.2 Cross-checking in other surveys 

In order to account for previous analysis, look for multi- 
epoch imaging data and additional colour constraints for 
the UCD candidates, we performed a cross-check with mul- 
tiple NIR and optical surveys, specifically the SuperCOS- 
MOS Science Archive (SSA), SDSS, Digitized Sky Survey 
(DSS), DENIS, WISE, the UKIRT Infrared Deep Sky Sur- 
vey (UKIDSS) and the SIMBAD database. We excluded 
nine of the selected UCD candidates, as they were classified 
as galaxies or fiare stars. Of the remaining 16 candidates, 
three were found to be known L dwarfs, another three have 
measured proper motions in the literature, and the remain- 
ing 10 do not have any proper motion measurements. This 
cross-check also allowed us to select the best epochs to cal- 
culate proper motions for the latter 10 UCD candidates. 
We used 2MASS as the first epoch and WISE as the sec- 
ond epoch when possible. Being the most recent NIR sur- 
vey, WISE allows us to have a larger baseline between the 
two epochs and thus more accurate proper motion measure- 
ments. If WISE images were not the most suitable ones, we 
then selected the second epoch from the other surveys. All 
images were chosen taking into account the largest base- 
lines possible and a reasonable number of appropriately po- 
sitioned bright stars that could be used as reference objects. 



2.2.3 Assessing common proper motion 

Proper motion calculations were made using standard IRAF 
routines. Measuring the positions of an average of 12 refer- 
ence stars, we have used the IRAF routine GEOMAP to 
derive the spatial transformation between the two epochs. 
Then the routine GEOXYTRAN applied these transforma- 
tions to the coordinates of the UCD candidate, and esti- 
mated how they changed between epochs. The uncertainties 
for this process are the combination of the rms scatter as- 
sociated with the coordinates transformation and the cen- 
troiding accuracies. Amongst our UCD candidates some had 
proper motions available in the literature. We compared our 
results with the published ones to check for consistency and 
then use the value with the smaller uncertainty. 

To verify if a pair has a common proper motion, we 
compared the proper motion in right ascension and declina- 
tion of each component, taking into account the uncertain- 
ties associated with each measurement. We consider that 
two objects have common proper motion if they match to 
within a < 2, a being the difference between their indi- 
vidual proper motions across the sky divided by the er- 
rors. In Table [5] we present the proper motions for all the 
candidates. Out of this sample, only two were found to 
have common proper motion with their primary star. The 
UCDs in these pairs are 2MASSJ133245. 30-^745944.1 and 
2MASSJ130541. 07+204639.4 (hereafter 2MASSJ1332-h7459 
and 2MASSJ1305-f2046 respectively). Information about 
these is shown in Tables [3] and 2] In one of the cases the UCD 
is an uncatal ogued L dwa r f, whi lst the other is classified as 
an L4 ± 2 bv lCruz et al. I (|2003t ). Our sam ple also inclu ded 
one known binary discovered by iFahertv et al. I l|2010l ). G 
62-33. 



2.3 Follow up observations 

The 2MASSJ1332-I-7459 UCD candidate was observed with 
the 3.58-m Telescopio Nazionale Galileo (TNG) in combi- 
nation with the Device Optimized for the LOw RESolu- 
tion (DOLORES or lrs for short) on 2011 May 13. lrs was 
equipped with a 2048 x 2048 E2V 4240 CCD and the lr-r 
grating with a resolving power of 7? ~ 700, centred at 7400 
A. A slit width of 1.0 arcsec was used with a seeing of ~ 1.1 
arcsec during the observations. This setup provided a wave- 
length coverage from 5200 A to 10400 A, with a dispersion 
of 2.7 A per pixel and a resolution of 3.8 A at the central 
wavelength. 

The spectra were reduced in the following way: the av- 
erage bias and fiat field correction was carried out using the 
FIGARO package from STARLINK using average bias and halo- 
gen frames, taken on the same night. The spectral reduction 
suite Pamela was used for the optimal extraction of the 
spectra. 

A neon and mercury arc lamp exposure at the same 
position as the targets allowed us to establish an accurate 
wavelength scale for each of the five spectra obtained. We 
used Molly, a 1-D astronomical spectral analysis package, 
to fit the arc frame with a fourth order polynomial, giving 
a rms of 0.17 A. The spectral resolution was measured by 
the full width at half-maximum (fwhm) of the arc lines, 
using an average of 10 lines located across the entire wave- 
length range. The individual spectra were flux calibrated us- 
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Table 2. Proper motion measurements for the L dwarf candidates 



Name of L dwarf candidates R.A. (J2000) Dec (J2000) pm^ ,4. pm^ec Possible primary Notes 



2MASSJ003142. 93-630818. 6 


00 31 42.94 


-63 08 18.6 


133 ± 8 


-88 ± 8 


HIP2540 


a 


2MASSJ071051. 38-492212.2 


07 10 51.38 


-49 22 12.4 


-14 ± 12 


121 ± 9 


HIP34739 


a 


2MASSJ074231. 27+180816.8 


07 42 31.27 


+18 08 16.8 


57 ± 31 


-78 ± 15 


HIP37622 


a 


2MASSJ100428.24-114648.9 


10 04 28.25 


-11 46 49.1 


-130 ± 15 


91 ± 15 


HIP49366 


a 


2MASSJ132434.95+545615.2 


13 24 34.94 


+54 56 15.4 


2 ± 12 


7 ± 10 


HR5055 


a 


2MASSJ134154.14-014553.1 


13 41 54.15 


-01 45 53.2 


-18 ± 5 


-34 ± 6 


HIP66886 


a 


2MASSJ134751. 79-104433.7 


13 47 51.79 


-10 44 33.7 


-432 ± 28 


-9 ± 18 


HIP67344 


a 


2MASSJ190536.28-370546.3 


19 05 36.29 


-37 05 46.3 


-0.2 ± 2 


-7 ± 6 


HR7227 


a 


2MASSJ073523. 28+315050.6 


07 35 23.18 


+31 50 50.6 


14 ± 10 


30 ± 8 


HD60179C 


a 



2MASSJ133245.30+745944.1 13 32 45.31 +74 59 44.2 -471 ± 28 39 ± 22 G255-34 b 

2MASSJ130541. 07+204639.4 13 05 41.07 +20 46 39.4 -23 ± 17 73 ± 27 GJ499AB b 



2MASSJ083204.51-012836.0 08 32 

2MASSJ132044.27+040904.5 13 20 

2MASSJ112149.24-131308.4 11 21 

2MASSJ022128.59-683140.0 02 21 

2MASS063447.73-582955.3 06 34 



04.51 


-01 28 36.1 


64 ± 13 


44.28 


+04 09 04.7 


-483 ± 19 


49.25 


-13 13 08.4 


-509 ± 10 


28.61 


-68 31 40.1 


46 ± 6 


47.73 


-58 29 55.3 


73 ± 9 



27 ± 15 GJ3504 c 

211 ± 17 HIP65121 d 

-81 ± 10 GJ3655 e 

-6 ± 17 HIPllOOl f 

-21 ± 9 HIP31300 g 



Notes: (a) New L dwarf candidates with new proper motion meas urements (b) new com mon proper motion pairs, here presented (c) 
prop er motion from Case well et al. ] ll2008h (d) pro per motion from 
bv lFahertv et al. I 1I2OI0I ) fe) proper motion from lSalim fc Gould 



Jameson et al. ] ll2008h and co mmon proper motion p air discovered 



motion from lRo'eser'et al. I 1I2OIOI ) 



. .. , pa 

1I2OO3I) ff ) prop er motion from lFahertv et al. I ll2009h (g) 



proper 



ing the nearby spectro-photometric standard flux star GRW 
+70d5824. The final spectrum is a variance weighted aver- 
age of the five individual spectra. 



3 BINARY CANDIDATES 
3.1 G255-34AB 

3.1.1 Properties of the L dwarf G255-34B 

To measure the UCDs spectral type, we have compared 
the optical spectra obtained at the TNG to standard 
template spectra f rom the L dwarf sequence defined in 
iKirkpatrick et al. I (ll999D . The comparison was done using 
four different subclasses, LO, LI, L2 and L3, using the follow- 
ing objects as standards: for LO, 2MASP J034532+254023; 
for LI, 2 MASSW J1439284+1929 15, for L2, Kelul (all 
three from IKirkpatrick et al.1ll999l '). and for L3, DENIS-P 
J1058. 7-1548 (|Delfosse et al. Ill997h . In order to define the 
best subclass for this dwarf, we used a m inimization and 
the g oodness-to-fit statistic, Gk, defined in ICushing et al~l 
l|2008f ). Whereas the minimization classifies this UCD as 
an L2 dwarf, the method described in lCushing et al."l (|2008l ) 
classifies it as an LI. The fit of the observed spectrum to the 
template spectra can be seen in Fig. [2] The fitting was done 
for the wavelength range 6300-9500 ^m. 

Secondly, we estimated three spec tral ratios, CrH-a, 
Rb-b/ TiO-b and Cs-a/VO-b, defined in IKirkpatrick et al. I 
l|l999l) for L dwarfs with spectral types < L5. The values 
for these ratios are 1.29, 1.08 and 0.91 respectively, usually 
associated with L2 dwarfs. 

Looking more closely at the spectrum of G255-34B, we 
see that despite showing TiO A8432, CrH A8611 and FeH 
A8692 lines all similar in strength, typical of LI dwarfs, the 
spectrum actuall y shows most of the spect ral key features of 
the subclass L2 (jKirkpatrick et al~lll999l ). As an example, 



we refer to the distinctly sloped part of the spectrum be- 
tween 7800 and 8000 A and the TiO lines A8432 and A7053 
that become weaker and disappear, respectively. 

Finally, we checked the two WISE bands to derive a 
spectral type in an independent way, using colour-spectral 
type pl ots. According to t he W l — W2 vs spectral type 
plot of IKirkpatrick et al~l (|201ll ) , the L dwarf, showing a 
Wl — W2 colour of 0.254, should have a spectral type be- 
tween M5 and L4. However, looking at the colours, J — W2 
and H — W2 vs spectral type plots, it can be seen that the 
UCD colours are similar to the ones shown by LO to L3 
dwarfs. These conclusions are consistent with the value ob- 
tained using our previous method. Taking all these spectral 
classification diagnostics into account, we here classify this 
UCD as an L2 dwarf. 

To determine the bolometric luminosity we have com- 
bined the optical spectra of the L dwarf with the photom- 
etry available. In the NIR part of the spectrum we used 
JHK photometry from 2MASS, extending to the mid in- 
frared with WISE ba ndpasses Wl and W2. We followed the 
method described in iMarocco et al. I l|2010l ) to estimate the 
bolometric luminosity. Firstly, we created synthetic spec- 
tra for those regions of the electromagnetic spectrum lack- 
ing observatio ns (NIR and mid-IR). We used the models of 
[Hubenv & Bu rrows 1 l|2007l ) since these cover a temperature 
range between 700 to 1900 K and have considered that the 
models are still valid for the typically larger temperature val- 
ues of early L dwarfs (up to 2250 K). We allowed for different 
values of log g, 4.5, 5.0 and 5.5, and assumed K^z = 10^, 10* 
and 10® cm~^s~^ (eddy diffusion coefficient). Secondly, we 
estimated the difference in flux between our models and the 
flux given by the available photometry, known as the scaling 
factor. Comparing this value with the scaling factor given by 
the distance assumed for the L dwarf and the radius range 
adopted (between 0.8 Kjup and 1.2 Rjup) we only consid- 
ered models for which the average scaling factor was within 
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Table 3. Properties of the system G255-34AB 




2MASSJ 1332+7459 



6000 7000 8000 9000 10000 

Wavelength (A) 

Figure 2. Optical spectrum of the 2MASSJ1332+7459 dwarf, 
with the standard spectra of an L2 dwarf, Kelu 1, and an L3 
dwarf, DENIS 1058-1548. The spectra has been displaced verti- 
cally in order to better compare them. The fit has been made for 
the wavelength range 6300-9500 lira. 



a certain interval. Thirdly, combining the synthetic spectra 
with the observed one, we calculated the bolometric flux, 
luminosity, and temperature range from the radius range. 
We then excluded the models in which the temperature was 
at least 100 K outside the constrained range. The final step 
takes into account only the remaining models and using the 
mean flux, we have estimated a final bolometric luminosity 
of 6.38 ± 0.67 X 10^^ erg s"^ or 1.67 ± 0.17 x 10"* Lq. 
Taking a radius range of 0.8-1.2 Rjup and the luminosity, we 
can derive the corresponding temperature for this L2 dwarf 
as Toff = 2080 ± 260 K. 

We have also derived the bolometric luminosity and T^s 
for the L dwarf following this procedure: taking into account 
the spectral type information derived before, we estimated 
th e K band bolometr ic corrections using the polynomial fit 
of iGolimowski et al. I ([2004) and hence the apparent bolo- 
metric magnitude of the object. Combining this with the 
distance, inferred from parallax measurements, we obtain 
the absolute bolometric magnitude and consequently, the 
bolometric luminosity. From the luminosity, T^s can be di- 
rectly calculated. The results are, for bolometric luminosity 
and Toff respectively, 7.14 x 10^^ erg s"^ or 1.87 x 10"" 
Lq, and Tes = 2106l^gQ K, and are in agreement with the 
previous values. For the Toff we again considered a radius 
range of 0.8-1.2 Rjup- 



3.1.2 System properties 

The components of the first system are the L2 dwarf dis - 
cussed in Section 3.1.1 and a K8 dwarf (jVvssotskv Ill956h . 
G255-34A, and a finder chart of these can be seen in Fig. |3l 
The spectral type we derived ea rlier for the L dwarf suggest s 
an Mj of 12.36 ± 0.11 based on lOupuv fc Liu et al. I l|2012h . 
This places it at a distance of 35.36 75 P'-i consistent with 
the measured parallax of the K8 primary, at 36.21 ± 1.50 



Parameter 


Value 


Name 


2MASSJ1332+7459 (G255-34B) 


R.A. (J2000) 


13 32 45.31 


Dec (J2000) 


-1-74 59 44.2 


J 


15.103 ± 0.041 


H 


14.116 ± 0.034 


K 


13.569 ± 0.032 


Wl 


13.170 ± 0.025 


W2 


12.916 ± 0.032 


Wl - W2 


0.254 


J - W2 


2.187 


H - W2 


1.200 


A^R. A. 


-471.296 ± 28.489 mas yr 




38.773 ± 22.113 mas yr 


SoTvDe 


L2 


lop" Ii 1 


-3.780 ± 0.045 


eit 


2080 ± 260 K 


Name 


G255-34A (HIP66074) 


R.A. (J2000) 


13 32 41.02 


Dec (J2000) 


-1-75 00 24.81 


B 


11.53 


V 


10.23 


R 


9.7 


I 


9.2 


J 


7.910 ± 0.024 


H 


7.302 ± 0.033 


K 


7.182 ± 0.016 


MR.A. 


-438.50 ± 1.33 mas yr "1 




49.82 ± 1.05 mas yr 


Parallax 


27.62 ± 1.14 mas 


Distance 


36.21 ± 1.50 pc 


SpType 


K8 


U 


-49.7 km s-l 


V 


-63.8 km s-i 


W 


-17.6 km s-i 


V total 


82.7 km s-l 


Teff 


4368 K 



pc (jvan Leeuwen IlioOTl l. G255-34A, a K8 dwarf, has a sepa- 
ration of 38.3 arcsec to the secondary, or 1364 AU. Initially 
believed to be part of a comm on proper motion pair with 
an other star by Luvten I ij 19791 '). it is not co nsider e d as s uch 



by^Cjriclas et al. 



197lll. as can be n oted in lWeisI (|l99ll ) 



ISchlaufinan fc Laughlin! (j2010l l have developed a way 
to calibrate M dwarf metallicity using photometry. Apply- 
ing their method and using V and Ks magnitudes, we find 
[Fe/H] to be —0.29 dex. However, metallicity uncertainties 
are significant for photometric constraints, ~ 0.2 dex. This is 
marginal evidence for G255-34A being slightly metal poor. 
Kinematics place this star outsid e the region defined as the 
Eggen box l|Eggen fc Iben |[l989h (and representative of the 
young disc population) and therefore suggest old disc mem- 
be rship and a lik ely age greater than 1.5 Gyr. According 
to lLeggett I ( | 1992|) this is broadly consistent with an [m/H] 
~ —0.5 and we employ the relationship for West's earliest 
spectral type range as reasonably representative for a K8 
dwarf to place a lower limit on the age of this system at 
0.8±0.6 Gyr. 
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Figure 3. A 4.5 arcmin X 4.5 arcmin finder chart for the system G255-34AB. On the left panel is the J band image and on the right 
one the W2 band. North is up and East is left. 



3.2 GJ499ABC 

3.2.1 Properties of the L dwarf G.J499C 

GJ499C is an L dwarf first discovered bv lCruz et al. I (|2003l ). 
The proper motion for this dwarf iias been estimated by 
some authors and tlie available values for this are shown 
in Table O Our own measurements are also presented and, 
within the error bars, are in g o od agreement w i th th e val- 
ues of IJameson et al. 1 l|2008l ). iFahertv et aTl (|2009l) and 
ISchmidt et al. I (|2010D . 

We have independently estimated the spectral type 
of the dwarf with the available spectrum from SDSS. We 
started by calculating the CrH-a, Rb-b/TiO-b and Cs- 
a/VO-b spectral ratios, with the additional colour-d ratio 
(since these are best suited to objects with spectral types > 
L5). The analysis of spectral ratios are all consistent with 
a spectral type of L5. We then compared the SDSS spectra 
with standard spectra templates, following the minimiza- 
tion and the factor method. We did this in order to assess 
if the fit to the spectral templates was consistent with the 
spectral type indicated by the spectral ratios alone. We find 
that minimization classifies the L dwarf as an L5, whereas 
the Gk factor suggests an L4 type. A final spectral type of 
L5 is a dopted, which is in agreement with the previous re- 
sult bv lCruz et aTl (|2003l '). that classifies it as an L4 ± 2. A 
finder chart for this L dwarf can be seen in Fig. |4l 

Following the method described in section 3.1.1, we 
have estimated the luminosity and T^a for GJ499C, using 
bolometric corrections and a radius range of 0.8-1.2 Rjup- 
The results are, for the luminosity and Tofi, 2.23 x 10^^ erg 
s-i or 5.83 X lO"'^ Lq, and 1574+"° 
what it is expected for L5 dwarfs 



and 1574j:}^^ K, in agreement with 



3.2.2 System properties 

The primary components of the syste m G499AB conta in a 
K5 with a closely separated M4 dwarf (iReid et al. I2OO4I). at 
a sepa ration of 0.9 arcsec. According to lDupuv fc Liu et al. I 
l|2012l ). the mean Mj value for a L5 dwarf is 13.56 ± 0.03, 
and with this we estimate a distance of 21.28lQ'2g PC to the 



GJ499 triple system. The parallajc meas urements from the 
Hipparcos catalogue (|van Leeuwen llioOTi ) show that the pri- 
mary is at a distance of 18.80 ± 0.61 pc, and is in agreement 
with the photometric distance derived for the UCD compan- 
ion. GJ499C has a on-sky separation of 8.6 arcmin from the 
primary, or 9708 AU. Prope r motion measurements come 
from the Ivan LeeuwerTI (|2007l ) catalog and are in agreement 
with other values publi s hed su ch as Zacharias ct al. (201^*), 
iKharchenko fc Roeser I (|20 0^ ) and IRoeser et al. I (|2008ll . In 
order to estimate the physical properties of this system, we 
follow the same approach as for G255-34AB. 

According to t he method described in 
ISchlaufman fc Laughlini (|2010l ) we can use the V and 
Ks magnitudes of 9.439 and 6.041 respectively (taken from 
the Hipparcos and 2MASS catalogues) for the K5 dwarf to 
estimate a [Fe/H] = —0.20 ± 0.2, suggestive of a slightly 
metal poor system. If we do the same calculations for the 
M4 companion, with V = 14.90 and Ks — 9.55, the result 
is [Fe/H] — —0.17 ± 0.2. Both results are consistent with a 

metal poor system. 

The velocity components published in iBobvl ev et al. I 

(2OO&) place the pair on the edge of the Eggen box, there- 
fore suggesting youth. Howe ver, kinematics offe rs a poor tool 
to place constraints on ages. IWest et al. I l|2008h used the ac- 
tivity lifetimes of M dwarfs to calculate ages, suggesting an 
age of 4.5l° Q Gyr for M4 dwarfs. We thus consider this as 
a lower limit to the age of the system. 



3.3 Possible contamination 

We performed a statistical analysis with the purpose of es- 
timating the probability that these two systems are, in fact, 
gravitationally bound and are not merely randomly aligned. 

Firstly we estimated the absolute magnitude in the J 
band for each UCD ± 1 spectral type. We did this using two 
different methods, first ly the mean 2MASS Mj as f unction of 
spectral type from the Dupuv fc Liu et al~l (|2012l 'l and also 
from the iMarocco et al. I (l201o | l rela tions between Mj and 
spectral type. Marocco et al. I j2010l ) uses the MKO photo- 
metric system, which we converted the 2MASS photometry 
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Figure 4. A 5 arcmin X 5 arcmin finder chart for the GJ499C dwarf. On the left panel is the J band image and on the right one the 
W2 band. North is up and East is loft. 



using the [Stephens fc Leggett! (|2004 ) relations. With Mj 
we estimate a minimum and maximum distance for each L 
dwarf. We then calculate a conical volume in the sky, with 
a radius that is defined by the separation between the UCD 
and its primary companion. For G255-34B the distances 24.5 
pc and 45.8 pc define a volume of 2.95xl0~^ pc^, whereas 
for the GJ499C, 12.6 pc and 27.0 pc define a volume of 0.116 
pc^ 

Secondly, we consider the numb er of stars within these 
volumes. The luminosity function l|Reid et al. 1 12007| ) de- 
rived from the 8 pc and 20 pc samples gives us a lower 
and upper limit for the density of stars in the sky (0.062 to 

0. 076 stars pc~^). Multiplying this number by the volumes 
around G255-34B, we expect to find between 1.830xl0~* 
to 2.242xl0~* MS stars that could appear as companions 
to the secondaries. For GJ499C we would expect between 
7.185x10"^ and 8.807x10"^ stars. 

Another important aspect to consider is that, even 
though stars could lie within the volume around the L 
dwarfs, only a small fraction could have common proper mo- 
tion. We therefore analyse the probability of finding Hippar- 
cos and Gliese stars that, by random chance, may have the 
same motion across the sky as our L dwarfs. In order to do 
so, we used the initial sample of MS stars with a maximum 
distance of 50 pc and considered firstly the ones lying inside 
a spherical volume with 45° radius. Such a large volume as- 
sures that we have enough stars in our sample to estimate 
the fraction of those that could appear as common proper 
motion to our secondaries. We have carefully excluded from 
this region a circular area where the common proper motion 
primary companion lies, to avoid counting it as one of the 
randomly aligned stars. Next, we selected only stars that 
could masquerade as common proper motion companions to 
the L dwarfs up to a 2a detection, which was used in our 
initial classification as a common proper motion candidate, 

1. e., within ± 120 mas yr~^ of the value of the secondary 
(valid for the two L dwarfs). 

Finally, taking into account the fraction of stars with 
similar proper motions as the secondary, and the sky den- 
sity, we find that for the G255-34B dwarf, an average of 



6.40 xlO~ MS stars could masquerade as a companion to 
the UCD by random chance, whereas for the GJ499C this 
value is of 1.36x10"'' MS stars. This statistical analysis al- 
lows us to say then, with a high degree of confidence, that 
both systems are genuine associations and not chance align- 
ments to a high degree of significance. 



4 DISCUSSION 
4.1 Binary fraction 

We have gathered information about all of the L dwarfs that 
are in binary or multiple systems, from the current litera- 
ture, finding 59 L dwarfs that form part of binary or multiple 
systems. Of these, 35 are in a system where the primary is 
another UCD. Two have white dwarfs as primaries and the 
remaining 22 have a MS star primary. For the purposes of 
this paper, we only present the latter 22, with information 
about these systems shown in Table[6l The distance has been 
estimated using parallax measurements where available or, 
otherwise, from photometric data. 

Considering the two new systems presented here, the 
sample has now 24 binaries, with six triples and two quadru- 
ples. This gives a fraction of 1:4 triple to binary systems, and 
1:12 quadruples to b inary systems. These values differ from 
the ones presented in lFahertv et al.1 (|2010l ). However, in the 
latter case the authors considered a wide UCD companion 
sample containing all UCD binary systems, and not only the 
ones with L dwarfs as secondaries. 

Our initial sample of L dwarf candidates from 2MASS 
was completed up to a magnitude of 16 in the J band. This 
is quite similar to the completeness limit of the Point Source 
catalogue of 2MASS. According to ISkrutskie et al. I (|2006h . 
2MASS is virtually complete (~ 99 per cent of the sky) to the 
10 a sensitivity limit for sources at or fainter than J = 15.8. 
Our sample is thus complete to this limit and we see that 
all known L dwarfs with J < 16 have turned up in our ini- 
tial search. Using the re lationship between absolute and ap- 
par ent magnitude (Marocco et al. l2010l . lKnapp et al. Iliooj 
and lLiu et al. II2006I ). for an L9 dwarf, we estimate that our 
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Table 6. Known binary systems containing one L dwarf secondary 



Name 



SpType 



SpTypc 



Distance (pc) sep(arcsec) projected sep(AU) Age (Gyr) 



Ref 
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(secondary) 


(primary) 










TT'Pi QO'7/1 A A 


r o 


T n 


QO /I Q _L n /I c 
oy.4o iL U.4e5 


Do 


z4oU 


z.l - I .Z 


/I 1 Q OC; "JQ 

4,lo,z0,oz,oo 




iVio 


T n 


iU. / o it U.lo 


9 cn 

Z.oU 


•^n 
oU 




/I 1 7 9c; 
4,1 / ,Z0 


ISJT TT997zl A 


iVl'rt 


T n 


oo.uu m z.uu 


9*^ n 

ZO.U 


4t(50 


A inn 

41:. - lU.U 


9*^ 
ZO 


A R Pir> A 


Ivl 


T 1 

IjI 


440. UD III l.^D 


0.0 


zoo 


n 09 ^'\ 9. 

U.UZ - lO.O 


'4:,O,Z0,Z0,0l 




IVvJ.O 


T 1 

1j1 


zl.Z/ it U.4o 


1 1 o 

11. y 


Z04 




A 1 9c; 
4, / ,Z0 




aivi4t.o 


T 1 _lT 1 


z / .4o it z. / U 


A A 
44 


1 A c;o 

140Z 




^^ ^(\ 9c; 
10,1D,Z0 






tip; 


1 /in nn -U nn 

i4U.UU it OU.UU 


n 7 

U. / 


yc5 


Q n 
o.u 


on OQ 
zu,zo 




iViZ 


J_jZ 


Qc: nn -U Q nn 
oo.uu it o.UU 


7Q n 
/ o.U 


loyo 


•3 A 
o - 4 


9Q 
ZO 




iVio 




1 A on -U 9 'rn 
14. yu it z. / u 
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Z40 
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U.UzO 


11 19 Qfi 
ll,lz,oD 


GJ 618. lA 


K7 


L2.5 


33.42 ± 3.00 


35 


1170 




4,13 


G 62-33A 


KO 


L3 


30.96 ± 0.82 


66 


2043 


4.4 


23,26 


G 200-28A 


G5 


L4 


45.66 ± 1.29 


570 


2600 


4.2 - 9.0 


23,26,27 


HD 49197A 


F5 


L4 


44.90 ± 1.21 


0.95 


44 


0.5 - 4.7 


4,9,26,27,28 


GJ 564A 


F9 


L4+L4 


18.17 ± 0.11 


2.64 


48 


0.05 - 12.4 


4,18,29,30 


HD 2057A 


F8 


L4+L4 


43.82 ± 1.69 


218 


9465 


1.6 - 5.3 


4,5,6,23,26,34,37 


Gl 417A 


GO+GO 


L4.5+L4.5 


21.93 ± 0.21 


90 


1953 


0.08 - 3.2 


14,15,27,29 


GJ lOOlA 


M3.5 


L4.5+L4.5 


13.01 ± 0.67 


18.6 


180 




1,2,3,25 


G 203-50A 


M4.5 


L5 


22.40 ± 1.90 


6.4 


135 




21 


LP 261-75A 


M4.5 


L6 


62.15 ± 28.58 


13 


450 


0.04 


10,12,25 


HD 203030A 


G8 


L7.5 


40.88 ± 1.24 


11.9 


487 


0.25 - 1.4 


4,22,25,26,28 


Gl 584A 


G0+G3 


L8 


17.86 ± 0.25 


194 


3465 


3.3 - 4.8 


19,26,34 


Gl 337A 


G8+K1 


L8+T0 


20.36 ± 0.22 


43 


875 


13.8 


13,24,26 


G 255-34A 


K8 


L2 


35.87 ± 149 


38.0 


1364 




38 


GJ 499AB 


K7+M4 


L5 


18.80 ± 0.61 


516 


9708 




38 



Ref ere nces: { 1) iGolimowski et al~ll2004l: (21 iMartm et al. Ill999l: (3) iHenrv et al. II2006I: (41 Undersoil fc Francis II2OIII: f5llLiu et al. I 



201C| : (6l lCruz et al. 112007 
199sl : (12)IShkolnik et al. iboOS 



(71 



Gizis et al II2OOII: (81 IChauvin et al. II2005I: (91 iMetchev fc Hillenbrand II2004I: (101 ?: (111 tRebolo et al. I 



(131 



Wilson et al. 



I I2OOII : 



14 



20031 : ( 1 61 ISeifahrt et al. ll20o'l 



Kirkpatrick et al. II2OOII: (15llBouvet al. .._ . , ._ . 

(17) 'Forveill e et al. II2OO4 I: (181 'Potter et al. "2OO 2I: (191 iMugrauer et al 11200 7^: (201 'Neuhauser et al. "200 ^1: (211 iRadigan et al. ll200Sl : 

(221 Met chev fc Hillenbrand 2pp6(@; (23 ) Fah ertv et al. II2OI0I: (241 iBur gasser et al. 2005^jg5rDupuv fc Liu et al. II2012I: (261 
[Casagrande et al. II2OIII: (271 iHolmberg et a l. 2009; (28) 'Writght et al. "2004; (29) Lafre niere et al. . 2007: (301 Lambert fc R eddv I 
I2OO4I : (311 iTetzlaff et aU 201ll (321 | Edvardsson et al. 1993 : (331 .Brvden et aL 2009; (34) ?; (351 ,Reid et al. ii200a : (36) 
IZapatero Osorio et ani201ol : (37) iReid et al.'TT2006bl : (38) this paper 



sample is complete up to a distance of 19 pc (for an LO this 
limit would be 94 pc). However, this completeness is only 
valid for the sky regions surveyed, since we have excluded 
reddened areas and avoided the galactic plane. 

In order to estimate the number of L dwarfs expected 
in the full sky, we took a spherical regi on of space with a 
radius of 19 pc, and multiplied it by the lCruz et al. I (|2007l ) 
space density estimate L dwarfs (3.8 ± 0.6 x 10~^ pc~^). 
We estimate that ^ 109 ± 17 L dwarfs should lie within the 
region studied. 

Based on Table [S] we find that six L dwarfs are part of 
binary systems in which the primary is a MS star, including 
the GJ 499ABC system presented here. Therefore, we find 
that the binary fraction for L dwarfs that have MS compan- 
ions is ~ 6 per cent. Considering now all the Hipparcos and 
Gliese stars, we find that 1787 were within our area and out 
to 19 pc. The binary fraction of stars with L dwarfs as wide 
companions is thus 0.33 per cent. The systems considered for 
this analysis have separations between 30 and 10000 AU. It 
is worth noting that the two systems with smaller on-sky 
and physical separations, G 239-25AB and GJ 564AB (with 
separations of 30 and 48 AU respectively), have been discov- 
ered with adaptive optics, since finding such faint objects so 
close to bright stars is rather difficult. 

If we now analyse the sample of known L dwarfs (based 



on Dwarf Archives), and do not take into account any dis- 
tance constraints, we see that the binary fraction of L dwarfs 
in any kind of system goes up to ~ 10 per cent. However, 
only 4 per cent of L dwarfs have MS stars as companions. 
We have not considered any possible unresolved binaries, 
and hence the values presented are just a lower estimate. 
The binary fractions here derived are not affected by the 
Malmquist bias, since both the number of single L dwarfs 
and the number of those found to be in multiple systems are 
affected in the same way. 



4.2 Unresolved binaries 

In order to identify any possible unresolved binaries in the 
sample presented in Table |6] we take advantage of the as- 
sociated high quality parallax measurements that can be 
inferred for the companion L dwarfs in known multiple sys- 
tems. We plotted the spectral type of the L dwarf com- 
ponents versus their Mj (Fig. [SJ. To estimate the abso- 
lute magnitude we used the primaries distances. We ex- 
cluded from this analysis two objects, GQ LupB (J15491209- 
3539039) and HD203030B (J21185897-h2613461) since they 
do not have 2MASS J band measurements. Also, we have 
plotted the G 124-62BC components as an unresolved bi- 
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Table 4. Properties of the system GJ 499ABC 



PcLrjimeter 


Value 


Name 


J1305+2046 (GJ 499C) 


R.A. (J2000) 


13 05 41.07 


Dec (J2000) 


+20 46 39.4 


J 


15.20 ± 0.053 


H 


14.04 ± 0.042 


K 


13.37 ± 0.039 


Wl 


12.540 ± 0.026 


W2 


12.153 ± 0.025 


Wl - W2 


0.387 


J - W2 


3.047 


H - W2 


1.887 


Mr. A. 


-42.61 ± 32.1 mas yr 




103.25 ± 15.4 mas yr 


SoTvDe 


L5 


loff Ij 1 Jjr\ 


-4.23 ± 0.05 


^ err 


1574 ± 170 K 


Name 


GJ 499AB {HIP63942) 


R.A. (J2000) 


13 06 15.43 


Dec (J2000) 


+20 43 44.40 


Vi 


9.40 


Ji 


6.89 ± 0.01 


Hi 


6.27 ± 0.01 


Ki 


6.05 ± 0.01 


V2 


14.90 


J2 


10.25 ± 0.10 


Hi 


9.75 ± 0.10 


K2 


9.55 ± 0.10 


MR. A. 


-49.17 ± 1.7 mas yr ~^ 


MDcc 


101.14 ± 1.1 mas yr 


Parallax 


53.18 ± 1.73 mas 


Distance 


18.80 ± 0.61 pc 


SpTypeA 


K5 


Sp TypeB 


M4 


T'eff 


4210 K 


U 


-9.0 km s-l 


V 


5.1 km s"-*^ 


w 


-1.4 km s~^ 


V total 


10.4 km s"i 


Notes: Vi, Ji, Hi and Ki refer to GJ 499A whereas V2, J2, H2 


and K2 refer to GJ 499B. 




Table 5. Proper motions for GJ 499C 


Ref. 


MR. A. (mas yr ^) /ioec. (mas yr ^) 


Zacharias et al. 2005 


-23.9 ± 9.0 -73.1 ± 9.0 


Jameson et al. 2008 


-23.3 ± 17.5 73.4 ± 26.9 


Fahertv et al. 2009 


-23.0 ± 17.0 73.0 ± 27.0 


Schmidt et al. 2010 


-59.1 ± 33.1 73.8 ± 22.6 


This paper 


-42.6 ± 32.1 103.3 ± 15.4 



nary, since there are no individual photometric measure- 
ments of the two LI components. The other five systems 
in which a MS star has two UCDs as secondaries are all 
resolved, and in Fig. [5] we plot the two components with 
different symbols, open circle and open triangle, and num- 
bered. The two new benchmark systems here presented are 
shown as filled circles. The filled squares are L dwarf com- 



panions that have not previously been considered as possible 
unresolved multiples. 

One interesting case is G 196-3B. The distance to this 
L2 dwarf has been e stimated by many different authors. 
IShkolnik et al. I (|2009| ) uses the companion, an M3 dwarf, to 
infer a photometric distance of 14.9 ± 2.7 pc to the system 
l|Reid et al. II2OO21 . 120071 ). iFahertv et al. I (|2009l '), on the other 
hand, uses 2 MAS S Mj vs spectral type relationships from 
ICruz et aL 1(1200311 to derive a different distance for this L2 
dwarf, 32.0 ± 2.0 pc. More recently. IZapatero Osorio et al. I 
(2010) discuss the system in detail giving minimum and 
maximum values for its distance of 15 and 51 pc, respec- 
tively. These two limits are given for two different scenarios, 
one if we consider the primary to be a single object and a 
field star, and secondly if we consider the secondary to actu- 
ally be a double object with ages of 3 Myr. Adopting a prob- 



able distance range of 15-30 pc (as IZapatero Osorio et ah 
2OIOI conclude and in accordance with the iFahertv et al. 



20091 and IShkolnik et al. I I2OO9I values), we estimate an Mj 
between 12.4 and 13.9. In[5]we have plotted the latter value, 
since this places the L2 dwarf further away from the poly- 
nomial fits. It can be seen that this particular dwarf shows 
re dder colours than expected fo r its spectral type, as noted 
in lZapatero Osorio et al.l l|2010l) (this is true even if we con - 
sider the dwarf to be an L3, as Zapatero Osorio et al. |[2010l ). 
The reason for this is still to be fully explained, although 
IZapatero Osorio et al. I (|2010| ) suggest the object might have 
a low gravity atmosphere with upper atmospheric layers or 
a warm dusty envelope. 

The L dwarf companion to AB PicA has been previ- 
ously reported as underluminous (|Fahertv et al. l2009t ). This 
is a young system, with ages smaller than 1 Gyr, and a low 
gravity dwarf. We note though that the unresolved binary G 
124-62BC, open square in Fig. O does not sit above the ex- 
pected position in the spectral type vs Mj plot. This could 
perhaps suggest that the actual single object sequence is 
slightly below the best-fitting sequences plotted, in which 
case AB PicB would not be underluminous, but making the 
position of HD16270B stand out even more as an overlumi- 
nous object. Just like HD16270B, LP 261-75B appears to 
be slightly overluminous for its spectral type. However, if 
we take into account uncertainties in the spectral type of ± 
1, we see that they both fit the sequence. 

Another interesting object is the system with HD2057A 
as primary, and two resolved L4 dw arfs as secondaries ( num- 
ber 2 in Fig. [5] Firstly mentioned in lReid et al~l (|2006bh . the 
two L4 dwarfs are late r confi rmed as companions to the F8 
star in iFahertv et al~l (|2010l ). As can be seen in the plot, 
the two symbols sit well above the average Mj for typical 
L4 dwarfs. In fact, one of the dwarfs h as Mj ^1.15 ab o ve the 
value given by the polyn omial fit of Marocco et al. 1 (|2010l ) 
and 1.01 above the one bv lFahertv et al. I (120121)7 This is not 
reported as a young system (ICa sagrande et ahlboill ). and 
therefore the overluminosity might be explained by multi- 
plicity. In such case, we could be looking at a quadruple 
system of L dwarfs. Only with further observations can we 
expect to truly explain this system. In the same way, the GJ 
564 system also appears to have one of the two L4 secon- 
daries above the sequence, and again binarity could be the 
explanation for its position in the plot. 
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Figure 5. Spectral type vs Mj for the L dwarf + MS star systems in Table |6] The filled squares are those systems with only one L 
dwarf, the filled circles are the two new systems presented in this paper. The systems in which there are two unresolved L dwarfs and a 
MS star are represented by two symbols, an open circl e showing the e a rlier o f the two UCDs and an open tr iangle showing the late r type 
UCD. The dashed line shows the polynomial fit from iF^hertv et al. 1 ||2012|'| and the solid one the fit from lMarocco et al. I 1 I2OIOI) . The 
crosses represent the mean 2MASS Mj for each spectral type jDuDuv~&: Liu et al. Il20l3 ). The numbers represent the following systems: 
(1) GJ 564ABC; (2) HD2057ABC; (3) Gl 417ABCD; (4) GJ lOOlABC; (5) Gl 337ABCD. 



5 CONCLUSIONS 

Using the all-sky surveys 2MASS and WISE, we have iden- 
tified two new systems with an L dwarf and a MS star as 
the primary components. Such benchmark systems allow us 
to better understand how the observed properties of UCDs 
are dependent on their ages and metallicities, and the two 
new systems here presented will help populate the growing 
sample of brown dwarfs with well constrained properties. 
Spectroscopy of the two primaries will be needed in the fu- 
ture to further estimate their ages and metallicities, thus 
fully characterizing the systems. 

We have also investigated the binary fraction for the 
specific type of binaries presented here, those in which the 
L dwarf has a FGK or M star as a wide companion. The 
number of such systems known to date is, however, small 
and any conclusions are dependent on the completeness of 
the sample, in this case, up to distances of 19 pc. The use of 
deeper surveys such as UKIDSS and VISTA will allow us to 
probe further distances and expand the number of such sys- 
tems. Considering that 2MASS is complete up to a J = 16, 
and taking into account the number of L dwarfs found, one 
can expect that by going three orders of magnitude fainter 
in the aforementioned surveys we can expect to increase the 
volume of sky sampled by up to 30 times. This means we 
can expect to find around 500 times more L dwarfs having 
a MS star as companion in an all sky sample. Ultimately, 
it is by studying a large sample of UCDs with abundances 
and age constraints from their companions that we can fur- 



ther understand and improve current atmoshperic models of 
these objects. 

Finally we show that some of the L dwarfs present in 
systems with a MS star could be unresolved binaries, due 
to their relative overluminosity for their spectral type. One 
interesting case is the HD2057 system, in which two resolved 
L4 dwarfs might actually be unresolved binaries themselves, 
thus making this a possible quadruple L dwarf system with 
an F8 star as the primary. 
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